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I. INTRODUCTION 


The need for lightweight, fuel efficient gas turbines 
that are capable of developing large amounts of thrust or 
power has motivated a continuing drive to obtain more ac- 
curate predictions of the flow through turbomachinery. 
Cascade testing of blade rows has, in the past, been a logi- 
cal and relatively inexpensive way to learn more about the 
phenomena involved in the flow through compressor and tur- 
bine stages. It is required more today in order to verify 
two dimensional and near-two dimensional analysis codes for 
flow through cascades. Such testing also provides two- 
dimensional blade element performance data which, in the 
absence of reliable analytical predictions, are required in 
the design of compressors and turbine stages. Reference l 
describes how cascade measurements are obtained using a cas- 
cade wind tunnel and then used in the design process. 

Before subsonic cascade wind tunnel data can be accepted 
as being valid, the flow conditions must meet three require- 
ments. These criteria are discussed in detail in Refs. l, 2 
and 3. First, any disturbance in the airflow should be 
caused by the test blades; that is, the inlet flow to the 
test section must be acceptably uniform. 

Secondly, the measured flow characteristics should, 


ideally, be independent of spanwise position along the test 


ie) 





blades. The flow, ideally, should be two dimensional. 

Duval [{Ref. 3] demonstrated that excellent flow conditions 
could be achieved in the Naval Postgraduate School Turbo- 
propulsion Laboratory (NPS/TPL) Subsonic Cascade Wind Tunnel 
using test blades with an aspect ratio of approximately two. 
The absence of suction along the walls results however, in 
some degree of streamline contraction which is measured in 
terms of an Axial Velocity Density Ratio (AVDR). 

The third requirement which must he satisfied is the 
periodicity of the inlet flow to the test section and of the 
outlet flow. Within one chord length of the leading edges 
of the test blades an upstream perturbation occurs as the 
streamlines adjust to negotiate the blade passages. Since 
the rectilinear cascade is simulating an infinite cascade of 
blades, the flow characteristics should be the same at cor- 
responding axial and blade-to-blade positions within each 
blade passage. This same condition should be true at any 
measurement plane downstream of the test blading. 

As described by Rose and Guttormsen [Ref. 4] several 
unique features were incorporated into the design of the 
NPS/TPL Cascade to ensure a two-dimensional and periodic 
flow at the test blading. Initial evaluations of the fa- 
cility were conducted and reported in Refs. 3, 4, 5 and 6. 
Work by Moebius [Ref. 6] involved modifications to the tun- 


nel plenum chamber which established satisfactory uniform 
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flow at the exit of the bellmouth contraction into the test 
section. 

In order to maintain an aspect ratio close to 2.0 ata 
solidity of 1.67 Cina [Ref. 7], following the work of Duval 
[Ref. 3], used a cascade configuration of 20 blades with 3 
inch spacing. Cina conducted a program of tests of DCA 
blading at five (5) different air incidence angles. ‘ith 
this cascade configuration, Cina found that the inlet flow 
to the test section was uniform in direction and of uniform 
static pressure, but with an imposed variation in velocity 
and stagnation pressure resulting from the wakes of inlet 
guide vanes. Although excellent periodicity was found over 
pairs of test blades, departure from strictly periodic con- 
ditions were detected from one blade passage to another. 
Cina explained this condition as being the result of the 
inlet guide wee wakes being separated at two inch intervals 
and entering a test section configured with a three inch 
blade spacing. Because of these flow conditions, Cina con- 
sidered his results to be preliminary. 

As a result of these findings the Cascade Wind Tunnel 
was modified so that inlet guide vanes were provided at one 
inch intervals. The object of the study reported herein 
was to obtain blade performance data on Cina's cascade with 
good periodic and uniform flow conditions. A necessary 
condition was to obtain agreement in the results for blade 


forces evaluated from surface pressures and from a momentum 


ony 





balance. A second objective was to compare measured blade 
surface Mach numbers with the results of code calculations. 

At the outset, it was first necessary to carry out an 
extensive testing program to verify the new inlet guide vane 
section and the effect the new spacing had on flow uni- 
formity and periodicity. It was found that the uniformity 
of dynamic pressure was improved with the inlet guide vanes 
Spaced at one inch intervals. Attempts were made to further 
improve the flow by the use of (various) wire screens placed 
downstream of the inlet guide vanes. These methods proved 
unsuccessful for the range of parameters tested and in fact 
aggravated the situation. 

Cina's testing of the Double Circular Arc blading was 
repeated without screens and with the Cascade Wind Tunnel 
configured with the modified inlet guide vane arrangement. 
Limited measurements were obtained before aero-mechanical 
problems with the new IGV arrangement, at the higher tunnel 
speeds, were encountered. 

The overall purpose of the testing program initiated by 
Cina was to obtain data with which to verify design optimi- 
zation computer codes developed by NASA. Towards this goal 
a fast, reliable computer analysis code (QSONIC) for cal- 
culating the flow field about a cascade of arbitrary 2-D 
airfoils was obtained from NASA. The code was adapted and 
modified to run on the Naval Postgraduate School's IBM 


3707 3033 computer. 
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The program QSONIC was developed by NASA to overcome the 
Mach number limitations of the earlier program TSONIC [Ref. 
8]. QSONIC is described in Ref. 9. Procedures for using 
the program QSONIC at the Naval Postgraduate School's com- 
puter facility are given in Appendix D. The procedures are 
documented for the case of the DCA blading in the NPS/TPL 
cascade wind tunnel. A program listing 1s included to docu- 
ment changes made to the code in order to adapt to the oper- 
ating system of the NPS computer. 

Preliminary results show that experimental measurements 


and code predictions are in very good agreement. 
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It. FACILITY DESCRIPTION AND MEASUREMENT APPROACH 


A. SUBSONIC CASCADE WIND TUNNEL 

The Naval Postgraduate School's Rectilinear Cascade 
Facility 1s shown in Fig. 1. A description of the facility 
as it was originally configured is given in Ref. 4. The 
test facility is an open cycle wind tunnel, designed for 
the purpose of testing cascades of axial-flow turbomachinery 
compressor or turbine blades. The unique design of the test 
section ensures that the airflow paths from the inlet guide 
vanes to all of the blades of the cascade test section are 
of equal length. This particular design was intended to 
eliminate the problems found in other cascade wind tunnels 
caused by having wall boundary layers of different thick- 
nesses entering the cascade at different points. 

As a result of the work reported in Ref. 1, two fine 
mesh screens were installed at the bellmouth entrance to 
improve flow stability. A follow-on study into the cascade 
performance was conducted by Bartocci and is reported in 
Ref. 5. As a result of Bartocci's findings, plenum turning 
vanes were installed to direct plenum inlet air towards the 
bellmouth entrance and to decrease the total pressure fluc- 
tuations present at the bellmouth entrance. Figure 2 shows 
the configuration of the plenum chamber as modified by Bar- 


tocci. Reference 6 describes work by Moebius that resulted 


20 





in further modification to the plenum chamber in which the 
original contraction was changed to two two-dimensional con- 
tractions in series. After this modification, acceptably 
small variations in velocity and flow angle were measured 

at the inlet guide vane station. Figure 3 shows the in- 
ternal arrangement of the plenum chamber as modified by 
Moebius and as it was configured for the work presented 
here. 

Using the plenum configuration shown in Fig. 3, Duval 
(Ref. 3] found that the wakes from the inlet guide vanes 
were not mixed out at the lower measuring plane of the test 
section but gave a well defined veriodic variation in the 
impact pressure. The peak-to-peak variation was +4% of 
dynamic pressure over two-inch periodic intervals. This 
condition was undesirable, but was tolerated while looking 
only to establish the values of parameters recuirec to 
achieve two-dimensionality and veriodicity. Since the 
inlet flow conditions were not uniform, mass averages were 
used to calculate properties at the inlet plane from probe 
measurements. 

In order to achieve a solidity of 1.67 anc aspect ratio 
of about 2, a blade spacing of 3 inches was required for the 
tests carried out by Cina [Ref. 7]. The tests showed unac- 
ceptable departures from blade-to-blade periodicity under 
conditions of high blade loading and the installation of 


additional guide vanes was recommended. The modification 
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to the inlet guide vane section of the tunnel resulting from 
Cina's findings is described in detail in Appendix A. 

In the present work, several tests were completed with 
the tunnel further modified by the introduction of wire 
screens between the inlet guide vanes and the lower plane of 
the test section. Appendix B Gescribes the screen material 
and the criteria used to select the particular screens used 


im thas study. 


B. INSTRUMENTATION 

The instrumentation used in the present study is that 
which 1s described in detail in Ref. 7. Twenty static pvres- 
Sure taps were located on the north and south side walls. 
The taps on the south wall were connected to a water mano- 
meter board so that the uniformity of the static pressure 
distribution of the inlet and outlet could be monitored 
visually. Additionally, one upstream tap on each wall and 
one downstream tap on each wall (near the centerline) were 
also connected to the Scanivalve so that these static pres- 
sures were recorded. 

Figure 4 shows the probe that was used for the upstream 
Survey (at the lower plane). The probe was a United Sensor 
Corporation DA 125 probe, serial number A847-1, calibrated 
earlier at various Mach numbers and yaw angles in a cali- 
Pratren facility. The United Sensor Corporation DC-125-24- 


F-22-CD probe, serial number A981-2 (Fig. 5), which was used 
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at the upper plane was similarly calibrated. The charac- 
teristics of the probes were approximated analytically to 
facilitate automatic data reduction procedures. The cali- 
bration and application procedures were those given by Duval 
(Ref. 3]. Appendix B of Ref. 3 describes both the upstream 
and downstream probes in detail. The mounting and traversing 


mechanisms are described in Ref. 7. 


C. REFERENCE MEASUREMENTS 

Plenum chamber (supply) pressure and temperatures, and 
atmospheric pressure were recorded on each data scan. 
Plenum pressure was also displayed on a water manometer 
board. The total temperature in the test cascade was as- 


sumed to be the same as the plenum chamber temperature. 


D. TEST BLADING 

The double circular arc test blading modeled the midspan 
section of the stator of the compressor stage reported in 
Ref. 10. Coordinates describing the profile of the blading 
are listed in Table D-2. The leading edge and trailing edge 
are shown in detail in Fig. 6. A photograph of the center- 
most blade is shown in Fig. 7. 

The three blades centrally located in the cascade were 
constructed with surface pressure taps along the midspan 
section as shown in Fig. 8. The centermost blade had 19 
ports on each of the vressure and suction surfaces and one 


tap at the leading edge. The two blades adjacent to the 
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center blade had 3 surface pressure taps located on each of 
the pressure and suction surfaces. The surface pressure tap 


locations for the centermost blade are given in Fig. 9. 


E. DATA ACQUISITION, REDUCTION AND ANALYSIS 

Data were recorded, reduced, and plotted using the modi- 
fied Hewlett Packard HP-3052A Data Acquisition System shown 
in Fig. 10. Reference 11 describes the system in detail. 
The system incorporated a HP-9845A desktop computer as a 
controller, with all components connected on the HP-98034A 
HP-IB Interface Bus. A NPS/TPL HG-78K Scanivalve Controller 
with two 48 port Scanivalves allowed the programmed acqui- 
Sition of probe and blade surface pressure measurements. 

The software used in the present study for acquisition, 
reduction and plotting of data were developed from software 
originally created by Duval and Cina. The programs are 
listed and described separately in Ref. 12. 

The uncertainties in the measurements are listed in 


Taele 1. 
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III. EXPERIMENTAL PROGRAM AND RESULTS 


=. “PROGRAM OF TESTS 

The test program was in three phases. First, in order 
to verify the new inlet guide vane assembly, tests were con- 
ducted with no blading in the test section and with the 
upper and lower endwalls set parallel at 35° (design condi- 
CLO) , 30° and 50° with respect to axial. 

Secondly, tests were made of the effect of wire gauze 
screen materials in reducing non-uniformities in the flow 
entering the test section. Appendix B describes the type 
of screens used and how they were installed. 

The last phase of the test program was a continuation 
of the work initiated by Cina. Table II lists the cascade 
configuration tested. One test was completed successfully 
before aero-mechanical problems were encountered and testing 


was halted until the causes were analyzed. 


B. 85ST PROCEDURES 
1. Cascade Adjustments 
In the first and second phases of testing, the same 
procedures were used to realign the cascade for each new 
configuration. The lower and upper end walls jere set to 
the desired flow angle and the inlet guide vanes were set 


so that their trailing edges were approximately aligned with 
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the end walls. The flow was started, and the desired inlet 
dynamic pressure was set. All tests were run at an average 
dimensionless inlet velocity (X) of about .13, corresponding 
to an inlet flow dynamic pressure of 18 inches water. Be- 
fore recording data, the water manometer board was checked 
to ensure that the distributions of wall static pressures 

at the inlet plane and outlet plane were acceptably uniform. 
If required, the inlet guide vanes were adjusted to obtain 
uniform static pressure to within +0.5 inches of water. 

In the third phase of testing, initially the pro- 
cedures used by Cina [Ref. 7] were followed, namely: the 
lower end walls were set to the desired inlet air angle and 
the upper end walls were set approximately to the expected 
exit air angle. The inlet guide vanes were set very approxi- 
mately and the cascade was turned on and set to an inlet dy- 
Mamic pressure of 18 inches water. The upper end walls and 
the inlet guide vanes were adjusted in turn to obtain wall 
static pressure distributions upstream and downstream which 
were acceptably uniform. Using this procedure however it 
was found on occasion that the inlet air angle sensed by the 
probe at the lower plane at mid-span could be 2 or 3 degrees 
different from the setting of the end walls. 

The following procedures was subsequently adopted. 
The lower end walls were set to the desired inlet air angle. 
The upper end walls were adjusted to be "wide open", to form 


a diverging passage in which, when the cascade was turned 
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on (to an inlet dynamic pressure of 18 inches water), the 
flow was completely separated. The inlet guide vanes were 
adjusted to obtain the required inlet air angle on the chan- 
nel center line over the center 24 inches in the blade-to- 
blade direction. The upper end walls were then moved 
individually towards the vertical until the lower vlane 
static pressure distribution was uniform and the upper plane 
Static pressure distribution was acceptably uniform at a 
value close to atmospheric pressure. No readjustment of the 
inlet guide vanes was made. 
2. Measurements 

Probe surveys were carried out in the blade-to-blade 
direction at midspan at the upper and lower planes. In the 
first and second phases of testing, data were taken over 
approximately 24 inches of the test section at intervals of 
0.25 inches. Also, in order to test the repeatability of 
measurements, repetitive samples were taken with the probe 
held fixed at midspan at the lower plane at the center, 10 
inches to the right and 10 inches to the left of center. 

During the third phase of testing, data were taken 


using the procedures established by Cina in Ref. 7. 


C. VERIFICATION OF INLET GUIDE VANE (IGV) ASSEMBLY 
The results of the first phase are presented (as shown 
img See Til) in Figs. 11 to 32. The results are arranged 


into groups. The first group (Figs. 1l to 14) are 
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measurements of tunnel conditions with the end walls set at 
35 degrees. Plots of conditions at the lower plane are fol- 
lowed by plots of conditions at the upper plane. 

Results for a wall angle of 30 degrees are given next 
(Fig. 15 to Fig. 18), followed by results for a wall angle 
me 20 degrees (Fig. 19 to Fig. 21). Data were taken over 24 
pmenes at the lower plane, and also at the upper plane at 
30°. At Sa) at the upper plane, only the center 12 inches 
were surveyed. 

The degree of repeatability of conditions in the wind 
tunnel from test to test (with no change in wall setting) is 
demonstrated by the results plotted in Figs. 22, 23 and 24. 

The last group of plots, Figs. 25 to 33, shows the de- 
gree of repeatability in the probe data from scan to scan. 
Data for these plots were obtained by holding the probe sta- 
tionary at midspan in three specific blade-to-blade loca- 
tions in turn and taking 50 repetitive scans of the channels 
normally recorded for survey profile data. The time inter- 


val for each scan was approximately 20 seconds. 


D. TESTING WITH WIRE GAUZE SCREENS 

The selection and installation of the wire gauze screens 
is described in Appendix B. 

The measurements obtained with the various screen con- 
figurations are given (as shown in Table IV) in Figs. 34 to 


45. The results are arranged in four groups. 
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Tae, farst group of plots (Fig. 34 to 40) give data ob- 
tained with the 16 mesh .0105 inch diameter wire screen in- 
stalled. Over a blade-to-blade distance from -1.0 to 7.0 
inches (Fig. 35) a peak-to-peak variation in velocity of 
about 1 percent was noted. This is slightly greater than 
the less than 1 percent (0.9 percent) variation noted over 
the same survey region without a screen installed (Fig. 12). 

The results shown plotted in Figs. 41 and 42 were ob- 
tained with two screens installed. One screen (16 mesh, 
-0105 inch diameter) was installed as discussed in Appendix 
B, while the second screen (2 mesh, .0400 inch diameter) was 
attached across the duct at the leading edges of the inlet 
guide vanes. 

The fourth qroup of plots (Figs. 43 to 46) show the re- 
sults for two different singhe screens. Probe survey data 
for these screens was taken only at the lower plane. The 
results shown plotted in Figs. 43 and 44 are data obtained 
with a 4 mesh, .041 inch diameter wire screen installed. 

The variation in velocity in the blade-to-blade direction 

was as much as t+1l.1 percent, neak-to-peak. The results 

shown in Figs. 45 and 46 are for a 5 mesh, .041 inch diameter 
wire screen. The variation in flow velocity was approxi- 


mately +1.5 percent, peak-to-peak. 
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E. PRELIMINARY TESTING OF DCA BLADES 

The results contained in Tables V to IX and Figs. 47 to 
60 are arranged in the following manner. 

The results shown plotted in Figs. 47 to 60 are divided 
into two separate groups. The first group (Figs. 47 to 57) 
contain results which exhibit the quality of the wind tunnel 
flow conditions. The second group (Figs. 58 to 60) shows 
the blade forces (and surface pressures) from survey data. 
In the first group of figures, results are presented first 
to examine the inlet flow uniformity (Figs. 47 and 48); 
second, to examine the outlet flow periodicity (Figs. 49 to 
53); and, finally, to examine outlet flow two dimensionality 
(ergs. S54 to 57). 


All points are shown connected with straight lines. 
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IV. DISCUSSION OF EXPERIMENTAL RESULTS 


A. EFFECT OF INLET GUIDE VANE (IGV) MODIFICATION 

The probe survey data shown in Figs. 11 and 12 were 
taken at the lower plane, with the end walls set at 35 de- 
grees. A turning angle of 35 degrees corresponded to the 
"design point" of the inlet guide vanes, when the airflow 
from the plenum chamber was at zero angle of incidence to 
the leading edge of the IGV's. These two figures show that 
the inlet plane total pressure at midspan in the blade-to- 
blade direction had a peak-to-peak periodic variation of 
about +2 percent and therefore about a +1 percent peak-to- 
peak variation in the velocity. Corresponding data from 
>probe surveys at the upper plane (Figs. 13 and 14) show that 
periodic variations in total pressure were reduced to about 
25% of the value at the lower plane by the mixing of the in- 
let guide vane wakes. 

Tt cam be seen in Figs. 15 to 21 that at “off-design” 
conditions for the IGV's (endwalls at 30° and 50°) there is 
a greater periodic variation in total pressure at the lower 
plane in the blade-to-blade direction than at the design 
point conditions. In Fig. 15 and Fig. 16, it can be seen 
that the periodic variations in total pressure are more pro- 
nounced with the flow from the vlenum at a negative incidence 


angle to the IGV's (endwalls at 0 ey. Except for the first 
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8 data points in Fig. 15, there is a well defined period of 
about 1 inch of travel. 

Figures 19 and 20 show the probe surveys conducted with 
endwalls at 50 degrees. At this off design condition the 
periodic variation in total pressure was considerably 
greater than the design point and the individual wakes from 
the inlet guide vanes were much less well defined. 

The repeatability of the survey was examined at wall 
angles of both 50° and 30°. Figs. 22 to 24 show that the 
non-uniformities in the flow conditions were reveated to 
(generally) better than 0.5 percent of total pressure. The 
guestion was then, to what accuracy could the individual data 
points be repeated in successive samples. This was examined 
at several probe positions and the results given in Figs. 25 


to 23 exolain the departures in Figs. 22 to 24. 


B. EFFECT OF WIRE GAUZE SCREENS 

All testing with screens was conducted with the end walls 
and inlet guide vanes set to yield a flow angle of 35 de- 
grees. The data obtained with screens installed were there- 
fore compared with the data obtained without screens, shown 
in Figs. 11-14. The effect of the pressure drop across the 
screen on the pressure coefficient plotted in Figs. 34-35 
should be noted. With the first screen installed the drop 
in total pressure from plenum to the probe in the lower plane 


was about 10 inches of water (plenum pressure minus total 
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pressure measured by the probe at the lower plane). With- 
out the screen (at design conditions), the pressure drop 
from the plenum to the probe at the lower plane was approxi- 


mately 2.0 inches of water. Since Q was defined as the 


ref 
difference between plenum pressure and lower wall static 
pressure, the value of QO -ef with the first screen installed 
was about 28 inches water and without the screen installed, 
about 20 inches water. In comparing the peak-to-peak varia- 
tion in Pl seen in Fig. 11 with that obtained with the first 
screen installed in Fig. 34, the difference in the values of 


Q must be considered. Examinations showed that the peak- 


ref 
to-peak variation in velocity remained at approximately 13% 
when the screen was installed. 

Figures 38 to 40 are plots of data obtained during a 
spanwise traverse of the probe at the lower plane. These 
figures show that the pressure drop through the screen and 
turning vanes was nearly uniform over approximately 8.0 
inches of the 10.0 inch span of the tunnel. 

In an attempt to generate upstream disturbances that 
might trigger early boundary layer transition on the IGV's 
and increase the rate of mixing of the wakes, a second screen 
was attached to the leading edge of the IGV's. The results 
in Figs. 41 and 42 showed that this was not the case and in 
fact the second screen increased the magnitude of the non- 


uniformities at the lower survey plane. Measurements made 


with a single 4 mesh (Figs. 43 and 44) and a single 5 mesh 
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screen (Figs. 45 and 46) of similar blockage showed that 
neither screen influenced the flow in a particularly favora- 
ble manner. The 4 mesh screen caused the peak-to-peak 
variation in velocity to be about +1.1 percent, while the 
> mesh screen caused the variation to be about +1.5 per- 
cent. This compared unfavorably with the variation ob- 
tained without any screen installed which was less than +1 
percent. 

It was therefore decided to proceed with measurements 


of the test cascade without using screens. 


C. PRELIMINARY TESTING OF DCA BLADES 
i. Imlet Uniformity 
The probe survey at the lower plane in Fig. 48 shows 
that the inlet plane total pressure at midspan varied in the 
blade to blade direction less than 0.5 inches of water, with 
no well-defined spatial period. This was an improvement in 
ee iniet conditions found by Cina [Ref. 7: Fig. 16]. That 
the spatial period was not well defined agreed with the 
findings presented earlier in this report. The wall static 
pressure distribution (Fig. 47) showed small variations 
(less than 0.5 ins. water peak-to-peak at the lower plane, 
.4 ins. water peak-to-peak at the upper plane). 
2. Two-Dimensionality 
Macmdata in Figs. 54 to 57 show that, at the down- 


stream plane, an area of (spanwise) nearly uniform conditions 
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existed near the centerline of the cascade. Reference 2 
points out that at higher loadings it is difficult to estab- 
lish a substantial spanwise area of uniform flow in the 
region near the suction side of the blade. This ai frvemiey 
is evident in the data shown in Figs. 54 and 55 which show 
that only about 20% of the spanwise distance is acceptably 
uniform. It is noted that Cina also found reduced areas of 
uniform flow at this incidence angle; however, 30~40% of the 
Spanwise distance was found to be acceptably uniform in his 
case. The difference could be the result of the reduced 
Spacing of the IGV's and its effect on the side wall boundary 
layers. 

Figure 58 shows results for inlet and outlet flow 
angles and blade force vectors derived in two ways as shown 
in Appendix B of Ref. 7. These two methods are first the 
applications of momentum conservation to probe survey data 
and second, the integration of surface pressures measured 
over the blade area. Reference 2 points out that for truly 
two dimensional flow the blade forces derived from the two 
methods should be the same. As shown in Fig. 58 the magni- 
tudes and the directions of the two vectors representing the 
blade forces are in reasonable agreement. It is noted how- 
ever that at this particular incidence angle Cina [Ref. 7] 
measured blade forces that were in total agreement in direc- 


tion but disagreed slightly in magnitude. The values of the 
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force magnitudes were about 1.5% lower than those measured 
by Cina. 
3. keeilodicity 

@5 Can be seen in Figs. 50 and 51 the total pressure 
and velocity qualitatively repeated fairly well over three 
central blade passages. Acceptably small quantitative dif- 
ferences are noted. There was also a small but measurable 
difference in the surface pressures on adjacent blades, as 
is evidenced in Fig. 49. 

4. Blade Performance 

Figures 59 and 60 are plots of the pressure and 
velocity distributions respectively over the centermost 
blade. These results compare favorably to those obtained 
by Cina for an incidence angle of Soe 

Table IX contains the blade performance parameters 
deduced from the probe survey data listed together with the 
data obtained for corresponding test parameters in Ref. 7. 
While differences in two sets of data are evident, the dif- 
ferences are not large. It is noted that the value of the 
loss coefficient was only 10% lower than was measured by 
Cina, but the AVDR was less than 2% different from unity 
rather than the 6.5% measured by Cina. Further measurements 
need to be made, particularly in the light of the following 


discussion, before stronger conclusions can be drawn. 
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2. Aero-Mechanical Problems Encountered 

Fifteen cascade tests were made while evaluating the 
new inlet guide vane assembly and testing the wire screens. 
All runs were made without test blades installed, with a 
plenum total pressure of about 20 inches of water and with 
the upper and lower end walls parallel. No difficulties 
were encountered in establishing the desired flow conditions 
or in using the inlet guide vanes to arrive at a satisfactory 
Gastribution of wall static pressures. 

The first time the Cascade Wind Tunnel was set up 
with test blades installed to take data at an air inlet angle 
of 26%..2°, the tunnel operated normally and the test was com- 
pleted. (The data from this test were subsequently found to 
be highly suspect and are not reported here.) During the 
next test, with an inlet air angle of 42.4, the start-up 
appeared normal and previously established procedures were 
used to arrive at a satisfactorily uniform wall static pres- 
Sure distribution. Tunnel operation appeared to be normal 
while taking data, but on shutdown a very noticeable high 
frequency vibration was encountered. Examination of the in- 
let guide vanes revealed that about 40% of the 60 blades 
were damaged. Damage included chips missing from the trail- 
ing edges, blades bent, cracks at the weld where the blade 
is joined to its support and indications that the suction 


side near the leading edge of one blade had been vibrating 
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against the pressure side near the leading edge of an ad- 
jacent blade. 

After the inlet guide vanes had been repaired and 
reinstalled extreme care was used at the beginning of the 
next test to adjust the inlet guide vanes, with two indi- 
viduals monitoring the movement of the adjustment mechanism. 
(IGV adjustment mechanism is described in detail in Appendix 
A.) <A lack of stiffness in the mechanism was suspected as 
having been a contributing factor to the failure. 

One successful test was completed at an inlet flow 
angle of 42.4° and these data were discussed above. 

With the next cascade configuration set, at an air 
inlet angle of ASG? when the IGV's were adjusted after 
starting up, high frequency vibrations were again experi- 
enced. The wind tunnel was shut down and no further testing 
was attempted at plenum total pressures as high as 20 inches 
of water gauge. 

The difficulty encountered with the IGV's is not 
fully understood, however the lack of stiffness present in 
the actuation of the two separate rows of vanes is suspected 
of having allowed the problem to occur. Certainly, the pos- 
sibility of an aerodynamic flutter condition being present 
(due to the misadjustment perhaps) can not be ignored. It 
was noticed after the initial failure that the lead screw 
which adjusts the IGV's could be turned but the blades 


mounted from only one side would be caused to rotate. This 
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could lead to the trailing edge of one blade contacting the 
tralling edge of an adjacent blade and effectively closing 
the blade passages. 

Also, the holes in which the cylindrical shanks of 
the IGV's were held were found not to be uniformly machined. 
As much as 0.1 inches of movement at the tip of some vanes 
was possible while others could barely move. (The most 
seriously damaged blades were found in or adjacent to the 
larger holes.) 

The tendency for the mechanism to "hang-up" on one 
Side would be greater as the vanes became more highly loaded. 
It 1s noted that the IGV problem was encountered first when 
going to increased incidence angles with the compressor test 
cascade installed. In setting a constant plenum total pres- 
sure of 20 inches of water gauge, the static pressure in- 
crease across the test blades to a constant atmospheric 
pressure at the downstream side implies that a progressively 
increasing dynamic pressure was being generated out of the 
turning vanes. This can be seen in Table IX, where QQ EO 


the test at 34 = 42.4 was 25 inches of water. 
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V. COMPUTATIONAL PROGRAM 


[= OSSCREPTION OF QSONIC 

The computational code, QSONIC, was developed by the 
staff at the NASA Lewis Research Center. This code is able 
to Calculate the blade-to-blade flow conditions in turbo- 
machinery blade rows assuming inviscid flow but including 
streamtube convergence and radius change in the throughflow 
direction. QSONIC is flexible enough to allow the input of 
the appropriate boundary conditions to calculate the flow 
through the test blading in the Subsonic Cascade Wind Tun- 
nel. The program uses a fully conservative solution of the 
full potential equation combined with the finite volume 
method on a body-fitted mesh. QSONIC uses an artificial 
density imposed in the transonic region, if such a region 
exists, to ensure stability and the capture of shock waves. 

The analysis used by QSONIC is a combination of transonic 
analysis methods to calculate the flow conditions in the vi- 
cinity of a cascade of airfoils. A conservative form of the 
full potential equation is discretized at every point of a 
body fitted periodic mesh and a mass balance is calculated 


through the finite volume surrounding the point. The volume 





Ithe helo and advice received from Charles Farrell at 
NASA Lewis R.C. in the process of adapting QSONIC to the 
NPS computer is gratefully acknowledged. 
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is corrected three dimensionally for any change in strean- 
tube thickness along a streamtube, if a quasi-3D solution is 
desired. Either elliptic or hyperbolic non-linear partial 
differential equations are used, depending on the local Mach 
number. 
The analysis used in developing QSONIC made the following 
assumptions: 
1) The airflow is inviscid and adiabatic. 
2) The airflow relative to the test blades is steady. 
3) Air is a perfect gas with constant specific heat. 
4) The airflow is isentropic and any discontinuities such 
as shocks are so weak that they may be approximated as 
HSCMEFOPLC JUMPS. 


5) There 1s no velocity component normal to the 
streamsurface. 


6) The airflow relative to a fixed reference frame (i.e. 
absolute velocity) 1s completely irrotational. 


Assumption 4 requires that the peak local relative Mach num- 
ber on a blade surface be 1.4 or less. The Mach numbers 
measured in test blades in the Subsonic Cascade Wind Tunnel 
would be well within this limit. However, this limitation 
would probably preclude the use of QSONIC for analysis of 
the flow field in the NPS transonic cascade wind tunnel. 
There are some combinations of blading geometry and flow 
conditions which cause unsatisfactory results to be gener- 
ated. For example, because of assumptions 1 and 6, sharp 
leading edges at high incidence angles (more that a few 


degrees) cause large velocity peaks in the blade surface as 
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the flow tries to turn from the stagnation point to the 
Suction surface. 

Reference 9 gives a detailed description of OQSONIC and 
the solution method used including the governing equations. 
Appendix D describes the operating procedures to use OSONIC 
on the Naval Postgraduate School's IBM 370/3033 computer. 
Appendix D also describes the input and output required as 


applicable to the Subsonic Cascade Wind Tunnel. 


B. APPLICATION TO THE TEST CASCADE 

Appendix D describes in detail the generation of the 
input required for QSONIC when applied to test blading in 
the Subsonic Cascade Wind Tunnel facility. In the present 
work one comparison of code calculations and measured data 
was made before teSting was stopped. The comparison was 
for an inlet flow angle (8,5) of 42.4°. 

Tables D.1 and D.3 show the input Gata generated. Table 
D.6 shows the flow solution output by QSOMNIC. Tne flow cal- 
culated on the blade surface, using a 15 by 97 mesh, was 
examined. Figure 61 is a plot of the calculated Mach number 
along the blade surface using two dimensional inputs. Fig- 
ure 62 is a plot of computed Mach number incorporating quasi- 
three dimensional effects. The method of incorporating 


quasi-three dimensional effects is explained in Appendix D. 
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C. COMPARISON OF CODE CALCULATIONS AND MEASURED DATA 

Table VII lists the data measured in the cascade wind 
tunnel. Col! C52 and Xe] are defined in Appendix C. The 
surface Mach number distribution measured on the center blade 
is shown plotted in Fig. 63. 

For comparison, the computed two dimensional, comouted 
gquasi-three dimensional, and the Mach number measured in the 
cascade wind tunnel are plotted together in Fig. 64. 

Excellent agreement between all three cases is seen. 

As would be expected, the greatest difference between mea- 
sured and calculated data is near the leading edge in the 


Suction side and at the trailing edge of the blade. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


Based on the results of the first part of the present 
Study, to evaluate the effects of the altered inlet cuide 


vane spacing on flow uniformity and periodicity, the follow- 


ing conclusions were drawn: 


1. With the inlet guide vanes operating at design, the 
peak-to-peak variation in velocity was about +£1%, and 


there was a well defined spatial period of about 1 
inen. 


2. Operating the cascade wind tunnel in a configuration 
that requires the inlet guide vanes to be set to other 
than zero incidence resulted in peak-to-peak varia- 
tions in velocity greater than +1% and a spatial period 
that was less well defined. 


The second part of the study, to evaluate the use of 
wire gauze screens to further reduce the non-uniformities in 
the flow field, led to the following conclusions: 

.l. A 16 mesh screen with a blockage factor of .69 had a 
Slightly aggravating effect in the variation in velo- 
city. The peak-to-peak variation in velocity with the 
screen installed was slightly greater than with no 
screen installed. This occurred at the expense of a 
pressure drop of about 10" of water across the screen. 

2. The use of screens with similar blockage but with 
larger mesh and larger diameter wires resulted in 
larger peak-to-peak variations in the velocity at the 
lower plane. 

The overall objective of the present study was to mea- 
Sure the performance of the DCA test blading. Because of 


aero-mechanical problems encountered with the inlet turning 


vanes the performance of the blades was obtained at one 
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incidence angle only. The following was concluded from the 


limited test program: 


18 


As a result of the reduced inlet guide vane spacing 
the variations in velocity and total pressure at the 
inlet plane were much less than those reported by Cina. 


Good periodicity was found from one blade passace to 
another. 


An acceptable region of spanwise uniformity (20-40% of 
blade span) was found at the downstream plane at the 
one test condition reported. However, this was less 
than was previously reported for the same incidence 
angle. 


The blade forces derived from the integration of sur- 
face pressure measurements and probe survey data were 
in close agreement in both magnitude and direction. 


The Mach number measured by surface pressure taps over 
the surface of the blade and Mach number calculated 
using the program QSONIC were in excellent agreement 
gualitatively and reasonable agreement quantitatively. 


The specific reasons for the aero-mechanical problems 
experienced with the inlet turning vanes have not been 
identified completely. 


Based on these conclusions and other observations, the fol- 


lowing recommendations are made: 


ie 


Use of the present inlet guide vane assembly and ad- 
justment mechanism for testing at inlet dynamic pres- 
sures higher than about 15 inches of water is unsafe. 
There are three possible solutions to this problem. 


a) Operate only within the dynamic pressure range of 
10-15 inches of water. 


b) Modify the new inlet guide vane assembly so that 
the vanes are supported at both ends on their axes 
of rotation. [Supporting the IGV's from both ends 
would prevent flapping vibrations of the (present- 
ly) cantilevered vanes. Such vibrations, when the 
vanes are supported alternately from opposite ends 
and the gaps are small compared to the chord, 
might lead to a potentially destructive flutter 
mode at particular flow velocities. ] 
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c) Replace the entire inlet turning vane section with 
one of entirely new design. 


The procedure should be adopted immediately of adjust- 
ing the vanes and walls of the cascade at lower values 
of the dynamic pressure before increasing the blower 
speed to the desired operating condition. 


More time needs to be spent to examine the flow field 
produced between the guide vanes and the test blades, 
and to establish the effects of the movement of the 
tall boards. The uniqueness of the flow field when the 
wall static pressure uniformity is used as a criterion 
of good inlet flow, needs to be examined by conducting 
repetitive tests at nominally similar test parameters. 
Only when the adjustment of the flow in the facility 
ane the quality of the flow itself is fully understood 
should the measured blade performance data be accepted 
as final. 


The upper electrical yaw adjustment mechanism should 
be replaced with a manual system to greatly decrease 
the time required to achieve probe pressure (angle) 
balancing. 


Develop the computer code necessary to take advantage 
of the plotting data created by QSONIC. 


Modify the data acquisition and reduction software for 


the HP 9845 so that real time plots of blade perfor- 
mance parameters can be displayed. 
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“ABLE I. MEASUREMENT UNCERTAINTY 


SS — Dr rsh SE SS see Spi cS > 


Item pescription Method Uncertainty 
—_—_— CC er et pi 


x Blade-to-Blade dimension Position ODay, 
x = Qin. West end Potentiometer 
x = 60 in. East end 








Z Spanwise dimension Position a 0) Eas Ly 9 
z= O ins North wall Potentiometer 
z= 10in. South wall on probe mount 
34 Inlet flow yaw angle Angle Potenti- +.2 deg. 


ometer on probe 
mount (hand 
adjustment) 





3, Outlet flow yaw angle Angle Potenti- eno deo. 
ometer on probe 
mount (motor 
driven adjustment) 


Plenum total pressure Statiestapein 01) ieee Og hs 8, 











3 len 4) 

z plenum chamber gauge 
V= 0 

P Static pressure at the Calibrated pneu- aol Dige 8 ere, 8, 
: test plane matic probe gauge 

E 2 Static pressure at Static tap on ae 007 aibele HO 
Ha ee in.,“y = -16.25 in., North wall gauge 

g- 0 an. 

E ™ Atmospheric pressure Absolute Strain ft . Seeins H,0 
a Gauge Transducer 

1 Pressure Scanivalve ot OL eeeer is H.0 


Transducer gauge 
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TABLE II. CASCADE CONFIGURATION FOR DCA BLADE TESTS 


Constant Parameters 


Number of Blades 20 

Spacing (Pitch) 3 inches 

Chord 5.01 inches 

Soduiidasty dae oY, 

Thickness 7.09 percent of chord 
Camber Angle 45.72 degrees 
Stagger Angle 14.72 degrees 


Variable Parameters 


34 42.4 degrees 


i 5.3 degrees 
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30 


50 


S10 


oi) 


30 


Tepe. iit. 


Survey Plane 


Lower 


Upper 


Lower 


Upper 


Lower 


Upper 


Lower 


Lower 


Lower 


SUMMARY OF MEASUREMENTS WITHOUT SCREENS 


Survey 


(24 


(24 


(24 


(24 


(24 


(Be 


(24 


(24 


Direction 


B-B 
inches) 


B-B 
inches) 
B-B 
inches) 
B-B 
inches) 
B-B 
inches) 
B-B 


inches) 


B-38 
inches) 


B-B 
inches) 


Fixed Probe 


@enen 


L. of ¢) 
(on €) 
Re of ¢) 


(anor 
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Pige 


SIL Is 3 


13 G4 


1D?-e 6 


BIR Bes: 


BES als B18 


25 


DIS OLE 


24 


Zo-2 7) 


28-30 


Bilas ye 


Nos. 


Purpose 


Flow Field 
Determination 


survey 
Repeatability 


Point 
Repeatability 





TABLE IV. SUMMARY OF MEASUREMENTS WITH SCREENS 


Aa Oo 
(By = 35°) 
SlliEy ey 
Screen Survey Plane Direction Fig. Nos. 
16 mesh Lower B-B 34 & 35 
.0105 wire 
Upper B-B 36 & 37 
Upper Spanwise 38-40 
i6 mesh Lower B-B 41 & 42 
.0105 wire 
+ 2 mesh 
aiWWead of IGV's 
4 mesh Lower B-B 43 & 44 
.041 wire 
5 mesh Lower B-B 45 & 46 
.041 wire 
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TABLE V. PROBE DATA, UPPER PLANE AT MIDSPAN (i = 5.3°) 


wAtTA FROM rile UWEREDS: Tie 
BLateS TO BLADE TRAVERSE MIDSPAH 


WEFERr PLAHE 


Sl i Or a Pita We ibar P=) (iba SE elieyar cet <a Ve 
BOSS TERE TCS TERETE S TERE SEECE RPE ER ee eee REESE HEE Re eK ETH ee See ECE OH 
1 a ee =]. 3 9656 « Zee t eS be oes 
e =¢. 54 =, oo i2o32 pala &, ~ led ec an! 
3 =e. ST =>, 38 242°6 2561 cose = cipal 
4 =Dinor 5 IS . 2S 2m? . 1 soe ml+3 
a =—S',,00 a2), 56 «dee. eee ~ Sem: me Bal ets 
G 4S =2 0 ~ar Sy este | ~Y3H2 a eal at 
i Dias = 1, ae . Die? Ree = 2316 « Sees 
3 -5.04 =1.. Sieerr coe 2.927 eats 
3 —4, oy =e Te Beet td sean 9x52 2.64398 
a3) -4,.58 ah, es sous 25938 ~¥9l oat 
11 7, a =i. =e « Joe fore ~y95%S Baris 
we —. 15 Se = ~ 9600 Pes) ~9954 otal 
13 =o. 3S -1. 35 woe? ~ 2636 Be a a oe 
14 =. a =. 35 9462 26382 e's oases 
iS = 39a =|. a 4053 . 2634 o lnae oa 
Le =. (39 =e. 33 - gfsde o one ~2263 Aboia te 
- =—6.. =a, ae 142a59 2) a 22455 edo st 
13 -3. 0% = Khe 24556 corse 16638 Pps) =o 
1? -2. 38 ie. at 2utor ~2r44 ~llws Ala SS) 
cas =—. ame =a, 9 a ara 2636 8962 ~5461 
os “a 9 =) ae Pe o\yl .2603 BS | cla »bdoe 
=2 2. 2 -1.F2 met tS ee oe 28926 ~o498 
26 sie), al -|.32 ~2814 fees ~WFL6 500° 
= =) ohare wk —-1. 32 ~731ld woe ou an) 6433 
Re, at Oe | —-1.468 wove Gi ~ 2606 ~Y 3701 To oue 
26 =. ae = 1h ie i oo Ome 2644 oo) 2 .8403 
ZV a) eee =e .974s 22526 ~asde ~b4r6 
23 ~ lee -1.46 .9orF »- 2661 .-lyud .ot+3y 
2+ = ans =| er "Oe 22°26 Pip ee re en -e 
30 =. ot -1{.46 » 4.794 wae 23 i gam S862 
oil -~.44 -1.4?7 . 4140 .29593 ~ 2698 ~3497 
32 — am - 187 8235 wZoac 22dde ae be be Te 
Bes = Om - 3.46 . 2811 .2rth ~1616 oak 
34 Biles. el ama . Ie pale ~l@dtek 2S 3c 1 
3 aoe =]. + 87ce “2602 8573 .o4ddd 
36 vow. =|. oro 2536 ~H365 2mo's 
37 .75 -1.33 . 5880 2591 HBT 5541 
38 a.. =|. a8 Pe | =) ee 213 alts ~6954 
3S Lame: =), Ss ~9902 oeal ~Hodt4 Pe = a 
43 ae =) 5 eae: BR ese sg FS Lee: 656] 
+1 | .oS =). oo wo ~2o8l W922 2 Seis 
$2 ie =l.33 4430 2 fae pene ner 
42 =. on -1. 34 4421 .267° wees o ott) 
+4 Je ~2. wl « Siete 2 ~ 2680 aS hee eto S 
4% Cis =e oe) ee ae ee ae ai ~YsVe2 2eti i 
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MABLE VI. PROEBE DATA, LOWER PLANE AT MIDSPAN (1 = 5.3%) 


TATA FRO FILE CRBREDS:T14 
BLADE Th SLADE TeaAverRte MI DSPAN 


LOWER PLANE 


oemt | 60Loc tim? Bet a We G1) bear Psiultar Pretiloar ope cha 

SHES FE CETTE STRESS EHHER EEC AETHER RHEER ERR REREEFERERE SSR SERS CH ERTS 
t = 4 He = HS oS er ee Pees: coe o heat e 
& =3. -42.42 0 PRIS -. 1444 239967 24rd 
3 - 3.9 “42.42 . 567 -. l4er7 296 22454 
4 es ~42.44 69614 =. to? ~ 1631 . 3458 
5 -2.00 -42,43 oe 4) -.1463 .O871 23524 
5 -1.Ge -42.43 o omar —, las 2.8875 ones 
v ~{. oe -42.43 . -. 1491 ~4GR'S oe 
= =." -42.42 9743 -. bos4 ~3921 22348 
3 8.0u ~42.43 eS =. Jone ees 235948 
1 . 96 ~-4$2,4¢2 ~ 3349 -,. 1S2e 2UgS4 . Sate) 
11 1.08 “$2.43 ~9814 -. io@a ~YdSSU Pate bell le 
Le t ."5e -42.43 . Me? 1 -,1628 O39 _ 2600 
i 2.0G -42.93 9937 =, [ou . 8390 2 Come 
14 Ce |) “42.9493 9742 -, 1622 2971 23950 
5 3.04 “42.41 . 9 7a -.1650 . 1859 ~33526 
is 3.59 =e. @- .98u¢ -. 1663 8924 AES) S) ot, 
LA 4.040 —~42.44 mee) oo =, 104 8968 acai 


Sv 





TABLE VII. CENTER BLADE DATE (i = 5.3°) 


Arlt mee lp gic: Mach mete | 


5 Sain i ie tinale it Te ina de ie ieia dh Rat Se ee ee ee oe a Salas Ei ie Ae ie Ei Ae i ie ale nie ie ne le i, ee eee 


preocesatMee SIDE CENTER SLADE 


eT ~Gut4 »c69> BS ae 5 lee PS ah 
~Ylea -00ls eoale? 4148 Le49 9824 
digo 230n5 area 2S A Ges eg meet ile 9951 
WETS Db) es soeor Bs 3 2126 9946 
2558 21S ~+5ri ~lser 2213 28953 
. Pairs 2583 4813 oles 2231 W395 
. PSs ~38452 247° RU elul? een . Lud 
i) 28576 24797 21120 seeds -H9a9 
2 3433 » 8603 47°64 21133 ~2243 ~O9 9c 
“132 28715 .4#892 . 1366 A Va 993 
2 $936 20736 .t37il «oer ccelg 2.9953 
ee . Oe? — 7 1 ~1146 e2idl -ooy 
other -0o73 249565 ates 221 PA ad (ee, 
0. LSe ~ Joel .453% - io els 0956 
r3s4 9487 ~ 4999 ~l14e3 fee . 3930 
2283 20411 .to/72 555 al ee . 1007 
inl ae Mere .4537 Urls oad ~1WSb 
eas 42-30 23245 ~913% 22354 old? 
22451 eer 3 235519 -.05'34 22445 . 1986 
73m ~ 89Ne oom -, 2452 2646 soleil 


I MEHTER BLADE 


xt 
iv) 


SUT Tord 


qn 


.0160 M@cS? -1. 3641 -3.2355 24973 Pivot 
.931% 8210 -.7756 =o eloay 24249 .L267 
ars .3329 -.49284 -1.6597 . 337s .176% 
.6853 .8553 -.3541 -1.3979 3875 leZe 
. bate .8712 -.3431 -1.37793 365% 21615 
.1956 .9976 -. 29% -1.2833 .3585 £1598 
. 2695 .117@ =, 2590 -1. 2684 _3Sz2 £1586 
63433 .1309 -, 2838 71.1243 .3453 21526 
14152 .1399 -,1540 -1.3338 .3377 21495 
£4938 .1432 - 3354 -, 9098 23276 ise? 
.5869 1412 -. 6554 -.3599 eee y. ices 
24a? . bee .8102 -.7351 Les J1ase 
.Tidts . 1269 .9631 » 7.6387 22s ol oes 
,Tres4 .1aZ1 .1456 -.4837 . 2286 J1298 
sae BESS .1911 -.4a59 .2385 en 
oe 9755 e326 -.3383 .27 38 es) 2 
ps2 24593 .2636 -,2748 ete s 1138? 
5 i , 3447 2882 -., 2365 .26a4 21178 
3530 B25 22748 -,2171 2613 .115!1 
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TABLE VIII. ADJACENT BLADES DATA (i = 5.3°) 


Aah Vem pl arses Nach tte | 


FST K SERRE SK SEER SORES H REE TEM RH HRT THREE EE RPS SHOE RTE MR RER ECR ETH RRETSE 


Pre-s'@E =s1DE LEFT BLADE 


7. ee aS Be .t014 =e 1 e442 BST) 

om eee . Biels -t019 9505 So 1o1 

3233 ee 24491 Wc so Pe) he . 162 
se tim SRE LEFT BLADE 

. 128 6°10 -,. 3433 =e se 360 ~lol 

41392 . 1406 —, Paes =1 8627 3 sol 1471 

wees 3 Wo9S Fis ~ 3294S Zr ry es 
PRESSURE SIDE FIGHT BLADE 

. Laos . S83 .to40 .890r Zar u - 161 

#1392 4716 4053 0933 2250 Lik 

sa203 O41 24597 8323 Bate, ~101 
SUCTIGH SIDE RIGHT BLADE 

213 U1 ~.3285 _ = sles 5506 Lew 
. 2 Poe 21409 =o Lee = leas > ~g30 4 143 
as25-3 W395 ae = om 22305 le4 
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4 fs 


a 


4 


CA fe iT 





TABLE IX. BLADE PERFORMANCE DATA 


Present Results From Ref. 7 

By 42.42 42.42 
a 54.5 5.3 
35 ge 0.4 
6 10.44 9.0 
D 05455 0.46 
W 0.087 0.041 
u) oe B 5 
[= 0.020 0.023 
20 cos By 
WJ cos B5 p 

oe (¥e0 ) 1.09 1.242 
AVDR L705 1.065 
CPomartc 0.413 0.351 
CoM -1.385 -1.380 
CoM ~0.669 ~0.566 

a ~1.476 

., es 30 
C =i 57 2 -0.645 
yB 
Q (in. H,0) Ze a 
x sie 512 
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Fig. 4. Lower Plane Survey Probe 
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Mig. 5. Upper Plane Survey Probe 
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Photograph of Instrumented Blade 
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Fig. 8. Instrumented Blade 
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Pig. ll. Probe survey Data at Midspan of Lower Plane 
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Fig. 12. Probe Survey Data at Midspan of Lower Plane 
Ema Walls at 35°, Points 51 to 100 
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Fig. 13. Probe Survey Data at Midspan of Upper Plane 
pagdemaiaks at 35°,;.Points 1 to 50 
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Pig. 14. Probe Survey Data at Midspan of Upper Plane 
Bad Walis at 35°, Bowmmts 51 to 100 
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Fig. 15. Probe Survey Data at Midspan of Lower Plane 
End Walls at 30°, Points 1 to 50 
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Fig. 16. Probe Survey Data at Midspan of Lower Plane 
End Walls at 30°, Points 51 to 100 


(Porenum ~ Pt) /2re¢ 


70 





Pip lenum)-Pl “Great 
eee ft Tt SS UPPER PLANE 3a O8G 





a 
| 
| 
| 
J Sa Jee) Ee Ree) Sn aed ESC, AE TE eats EOL ee O_O Cn ekenrt eee [oe 
od hog oo ul ag ny 
~ — 1 | | | 
! t 
iene = 


Fig. 17. Probe Survey Data at Midspan Of Upper Plane 
End Walls at 30°, Points 1 to 50 
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Fig. 18. Probe Survey Data at Midspan of Upper Plane 
End Walls at 30°, Points 51 to 100 
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Fig. 19. Probe Survey Data at Midspan of Lower Plane 
End Walls at 50°, Points 1 to 50 
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Fig. 20. Probe Survey Data at Midspan of Lower Plane 
Bad Walls at 50°, Points 51 to 100 
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Fig. 21. Probe Survey Data at Midspan of Upper Plane 
End Walls at 50° 
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Fig. 22. Probe Survey Data at Midspan at Lower Plane 
End Walls at 50°, Two Runs, Points 1 to 50 
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Fig. 23. Probe Survey Data at Midspan at Lower Plane 
Emad Walls at 50°, Two Runs, Points 51 to 100 
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Fig. 24. Probe Survey Data at Midspan at Lower Plane 
End Walls at 30°, Two Runs, Points 1 to 50 
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Fig. 28. Repetitive Samples with Fixed Probe Position 
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Fig. 29. Revetitive Samples with Fixed Probe Position 
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Fig. 30. Repetitive Samples with Fixed Probe Position 
(On Centerline at Midspan, End Walls at 30°, 
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Fig. 33. Repetitive Samples with Fixed Probe Position 
(10° Raght of CTR Midspan, End Walls at 30°, 
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Fig. 34. Probe Survey Data at Midspan, Lower Plane 
be Mesh Sereen, Walls at 35°, Points 1 to 50 
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Fig. 35. Probe Survey Data at Midspan, Lower Plane 
16° Mesh Sereen, Walls at 35°, Points 51 to 100 
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Fig. 36. Probe Survey Data at Midspan, Upper Plane 
16 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 37. Probe Survey Data at Midspan, Upper Plane 
16 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 41. Probe Survey Data at Midspan, Lower Plane 
16 Mesh and 2 Mesh, Walls at 35°, Points 1 to 50 
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Fig. 42. Probe Survey Data at Midspan, Lower Plane 
16 Mesh and 2 Mesh, Walls at 35°, Points 51 to 100 
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Fig. 43. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 44. Probe Survey Data at Midspan, Lower Plane 
4 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 45. Probe Survey Data at Midspan, Lower Plane 
9 Mesh Screen, Walls at 35°, Points 1 to 50 
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Fig. 46. Probe Survey Data at Midspan, Lower Plane 
5 Mesh Screen, Walls at 35°, Points 51 to 100 
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Fig. 48. Probe Survey Data at Upstream Midspan 
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Fig. 49. Blade Surface Pressure Distribution on 
Three Centermost Blades (i = 5.3) 
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Pig. 51. Probe Survey Data at Midspan 
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Fig. 53. Probe Survey Data at Midspan 
(1 = 5.3, Outlet Angle, Upper Plane) 
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Fig. 54. Spanwise Probe Data Surveyed 1 in. from 
suetion Side of Centermost Blade 
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Fig. 55. Spanwise Probe Data Surveyed 1 in. from 
Suction Side of Centermost Blade 
eos .5, </X, Upper Plane) 
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Fig. 56. Spanwise Probe Data Surveyed 1 in. from 
Pressure Side of Centermost Blade 
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Fig. 57. Spanwise Probe Data Surveyed 1 in. from 
Buessure Side Of Centermost Blade 
(1 = 5.3, X/X, Upper Plane) 
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Fig. 58. Resultant Blade Force Vectors by Momentum 
Balance (-~--—--— ) and from Surface Pressure 
Integration (————) i= 55S 
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Fig. 59. Measured Blade Surface Pressure Distribution 
(1 = 5.3, * = Pressure Side, 
+ = Suction Side) 
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Fig. 63. Measured Blade Surface Mach Number Distribution 
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APPENDIX A 


MODEPICATION TO THS INLET GUIDE VANE SECTIOM 
OF THE SUBSONIC CASCADE WIND TUNNEL 


As discussed in Section I, Cina discovered during his 
test program that while the inlet flow to the test section 
was uniform in direction and uniform in wall static pres- 
sure, it contained a variation in velocity and stagnation 
pressure resulting from the wakes of the IGV's. Because the 
inlet guide vanes were spaced at intervals of two inches and 
the test cascade blades spaced at three inches, departures 
from strictly periodic conditions were detected from one 
test blade passage to another. 

To alleviate this problem, the inlet guide vane arrange- 
ment was modified so that the guide vanes were placed at l 
inch intervals. In order to preserve the option of reverting 
to a two inch IGV arrangement and because it was not possi- 
ble to machine the south wall to hold additional blades, a 
separate structure was placed between the bell mouth contrac- 
tion and the walls of the cascade. By mounting the IGV's 
in a separate unit which remained fixed once installed, 
hardware adjustments between tests associated with a change 
of end wall angle were greatly simplified. 

The new inlet guide vane assembly was constructed using 


two lengths of 10 inch steel channel as shown in Fig. A.l. 





One set of guide vanes was mounted on the south side of the 
unit at 2 inch spacings. A second set of vanes was mounted 
at 2 inch spacings on the north side of the unit. When the 
unit was assembled the guide vanes mounted from alternate 
sides meshed, resulting in an inlet guide vane spacing of 

1 inch. The unit was provided with a single hand crank at 
the east end so that the vanes would be adjusted in unison. 
Once installed the complete structure could be left in place 
when the cascade north wall was removed to adjust air inlet 
angle. The one inch vane spacing ensured that periodicity 
at the test section would result for any test blade spacing 
which was a multiple of 1 inch. Equally important, the 
wakes remaining at the inlet to the test cascade would be 
greatly reduced as a result of closer spacing. 

Figure A.1 shows the details of the inlet guide vane 
unit, while Fig. A.2 shows the assembly in relation to the 
bellmouth contraction and the side walls. Figure A.3 shows 
the mechanism to adjust the inlet guide vanes. Figure A.4 
shows a view of the IGV assembly from the north side. A 
view of the Cascade Wind Tunnel partially assembled (north 


Side wall off) is shown in Fig. A.5. 


0 








Fig. A.1. Inlet Guide Vane Assembly 
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View of the IGV Adjustment Mechanism 
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View of the Subsonic Cascade Wind Tunnel (North Wall Removed) 
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SELECTION AND INSTALLATION OF SCREEN MATERIAL 


Pankhurst and Holder [Ref. 17] show that the turbulence 
of an airstream can be increased by placing a coarse mesh 
across the flow upstream of the test section. One of the 
most effective methods which is used for the reduction of 
turbulence and non-uniformities also consists of placing a 
mesh screen across the tunnel. Screens used for this pur- 
pose are of a finer mesh and are placed at a greater dis- 
tance from the test section, and normally in the low-speed 
region upstream of the contraction in a conventional sub- 
sonic wind tunnel. By using such a screen the large scale 
eddies are removed at the expense of the introduction of a 
greater number of smaller eddies which decay rapidly. 

McEligot [Ref. 16] investigated the possibilities of 
reducing non-uniformities in the test section of the sub- 
sonic cascade wind tunnel. His investigation and recom- 
mendations were completed while the inlet guide vanes were 
still at 2 inch spacings. McEligot concluded that some 
modification was necessary to achieve one percent uniformity 
for the mean velocity at the test cascade inlet plane and 
Suggested several options. One of the options suggested was 
placing the turning vanes (inlet guide vanes) at a closer 


pitch. As explained in Appendix A, the pitch of the inlet 
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guide vanes was reduced from 2 inches to l inch. This new 
inlet guide vane arrangement did result in a one percent 
uniformity for the mean velocity at the test cascade inlet 
plane. 

The other approach suggested was the use of wire gauze 
screens. McEligot showed that the velocity distribution ap- 
peared to be largely dependent on a pressure drop coefficient 


K, defined by the equation 


Py [8P5 


so V7 


where Py and P5 are the pressures upstream and downstream of 
the screen respectively. This pressure drop coefficient de- 
pends mainly on the blockage coefficient 8 defined by the 


equation 
B= (1 - a/2)? 


where d is the diameter of wire used in the screen and 2 is 
the distance between the wires. This blockage coefficient 
is commonly referred to as "percent open area" in catalogues 
of industrial wire cloth and woven wire screens. 

For the velocities and flow angles used in the cascade 
wind tunnel, McEligot recommended using a wire gauze screen 
with a resistance coefficient, K, of 2.2 and a blockage 
coefficient, 8, of 0.47. However, since the new inlet guide 


vane arrangement resulted in a one percent uniformity for 
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the mean velocity and the pressure drop across a screen with 
a blockage coefficient of 0.47 was expected to be higher 
than could be tolerated for the desired test conditions, 
screens with a slightly higher blockage coefficient (higher 
percent of open area) were selected to be tested. 


The screens tested were of the following configurations: 


MESH WIRE DIAMETER (inches) BLOCKAGE COEFFICIENT 
+ .0410 ~ 69Is9 
5 .0410 S15) 870) 
16 G20 5 sOu22 


Until the effectiveness of wire gauze screen in reducing 
non-uniformities in this facility was proven, a temporary 
installation of the screen material was considered adequate 
for testing purposes. The test screen was installed in the 
cascade wind tunnel by placing it between the inlet guide 
vane assembly and the north and south end walls. This ar- 
rangement placed the screen 7.25 inches downstream of the 
inlet guide vanes and19.3 inches upstream of the lower test 
plane. Figure B.1l shows the installation of the wire gauze 


screen. 
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APP ENDEX C 


CASCADE PERFORMANCE PARAMETERS 


(em F. S. Cina; reproduced with minor changes from Ref. 7) 


The performance of a cascade 1s specified in terms of 
the deviation angle (5) and the loss coefficient (w) for 
given inlet conditions. In Ref. 1 the loss coefficient is 
shown to correlate in terms of the Diffusion Factor (D). In 
the present work, the performance parameters were calculated 
using the following expressions: 


1. Loses Coefficient (w) 


o = ——— (1) 


where the mass averaged pressure coefficients in Eq. (1) are 
defined in Appendix C of Ref. 7. It is shown in Appendix 

C of Ref. 7 that the effect of time dependent supply con- 
ditions are removed and the effect of AVDR is included in 
the use of Eq. (1). 


2. Diffusion Factor (D) 


WwW. * Dow (2) 


aR ei 0) 





> where VY is the local velocity, V 


3. Pressure Rise 


C = bos (3) 
Pstatic oF 


4. Blade Surface Pressure Coefficients 


GC => (4) 


Q 
II 


(5) 
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velocity and T, is the stagnation temperature. 
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APPENDIX D 


INSTRUCTIONS FOR PREPARING INPUT AND OPERATING 
QSONIC USING A RECTILINEAR CASCADE CONFIGURATION 


B.1 BACKGROUND INFORMATION 

QSONIC has the capability to calculate an axial, mixed 
or radial flow field and the test cascade can be rotating or 
stationary. The geometry of the streamsurface can be a 2D 
planar cascade or axisymmetric with varying channel thick- 
ness and radial position. The capabilities of QSONIC, be- 
yond those of previous cascade analysis methods (such as 
described in Ref. 8) include the ability to calculate 
through weak shocks with a peak relative Mach number less 
than 1.4, and completely around both leading and trailing 
edge regions of a blade profile. The blade shapes in the 
leading and trailing edge regions are not restricted to 
Circular arcs. Detailed instructions for preparing input 
for a configuration other than an axial flow, stationary and 
rectilinear cascade may be found in Ref. 9. What follows 
are instructions for preparing the input applicable to the 
Rectilinear Cascade facility and running QSONIC on the Naval 
Postgraduate School's computer and associated operating sys- 
tem. It is assumed that the reader has a working knowledge 
of the NPS computer operating system and is familiar with 


Refs. 13 and 14. 
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QSONIC operates in two parts. The first part gener- 
ates a body (blade) centered mesh (geometry generation). 

The second part actually solves for the flow conditions at 
points in the mesh (flow solution). The data necessary to 
generate a mesh consists of a two-dimensional description of 
the blade shape. This is in the form of pairs of (X,Y) 
points on the surface, together with parameters that de- 
scribe the cascade layout, such as chord and stagger angle. 
Additionally, parameters describing the density of mesh 
lines complete the input for the geometry generation. 

For flow field calculations, the upstream flow condi- 
tions, convergence criteria and a schedule of meshes to be 
used should be input. If quasi-three dimensional effects 
are to be considered, a data file containing a description 
of the streamchannel's radial thickness and position as a 
function of distance along the stream surface is needed. 

For the case described herein, this was input by assuming a 
linear reduction in streamchannel thickness using a factor 
of 1/AVDR. This gave excellent results. (The output of 
another NASA code, Meridl [{Ref. 15], can be used to input 
data to QSONIC for compressor flow field calculations. This 
program has recently become operational on the NPS computer.) 

The output of QSONIC consists of listings which con- 
tain an echo print of the input data, generated mesh coordi- 
nates on the blade surface, progress reports on the flow 


convergence and a list of the final velocities, pressures 
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and densities on the blade surface for each grid that was 


gnclu@ed in the schedule of solution meshes. 


Deve INPUT DESCRIPTION 

The input for QSONIC falls into the following 
categories: 

Logical and case control parameters (NAMELIST PARAMS) 
Bulk data input: 

For geometry generation runs (NAMELIST INSTUFF) 

For flow solution runs: 

Mesh point storage files 

Streamchannel data file (for quasi-3D) 

The following is a description of the logical and case 
control parameters for the Rectilinear Cascade. Except for 
TITLE, the format for all these variables is in namelist 
form. The namelist is PARAMS. This information is taken 
from Ref. 9 and adapted to the Cascade Wind Tunnel. Name- 
list variables can be entered in any order. If a default 
value is listed, it is not necessary to enter that particu- 
lar variable. If the default value is listed as none, then 
a value for that parameter must be input. 

The following parameters apply to both the mesh gen- 
eration and the flow solution runs, but the values need not 


be the same. 








Variable 


Name Type 
TETLE Alphameric 


NOFLOW Logical 


MS Integer 


NOZES Integer 


Default 


Value 


None 


False 


None 


None 


Description 


This is a one-line name for 
thes case bermg run. TITLE 
Must appear on the first line 
of the data file that will be 
referred to as NAMELIST DATA. 


NOFLOW = .TRUE. if a run is 

to stop after generating a 
mesh, such as the mesh genera- 
then run. 


NOFLOW = .FALSE. for the flow 
Solution wun. 


MS iS an array (max. dimension 
10) of values (max. value = 
2a) fer “eke numiser of grid 
lines in the mesh that will 
enclose the blade. 25 1S a 
satisfactory number for blades 
of solidity near unity. As 
solidity increases the maxi- 
mum value in MS should de- 
crease. 


NOZES is an array (max. 
dimension 10) of values (max. 
value = 49) for the number of 
grid lines in the mesh ra- 
diating from the blade on one 
surface. If a value of NOZES 
1s greater than zero, a mesh 
with that many lines will be 
developed and stored in the 
file MESHGEN DATA. If the 
value is negative, QSONIC 
assumes that this MESH al- 
ready exists in the file 
MESHGEN DATA and will be read 
in. For Geometry Generation 
runs NOZES > O and for flow 
Solucion curs NOZES < 0. For 
electrostatic analog grid 
generator, NOZES must be odd. 
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BETA] 


BETA2 


GAMMA 


TOES 


Real 


Real 


Real 


Real 


None 


None 


ena 


None 


BETA] is the flow angle at the 
upstream boundary, in degrees. 
BETA1 is measured from the aero- 
dynamic chordline to the direc- 
tion of flow; clockwise is 
negative. This can be obtained 
from the output of the cascade 
tunnel data reduction program 
'CX4431'. It is listed as "“in- 
let air angle." Because of a 
difference in definition, it 

1s necessary to use the usual 
inlet air angle minus the stag- 
ger angle for BETAIL. 


Example: The output of CxX4431 
lists an inlet air angle of 
42.429°. The cascade is con- 
figured with a stagger angle of 
4. 27. 


a2 .429° 
- 14.27 


28.16 = BETAIL 


BETA2 1s the flow angle at the 
downstream boundary in degrees. 
It 1s measured from the aero- 
dynamic chordline to direction 
of flow. Clockwise is negative. 
This can be obtained from the 
output of the cascade tunnel 
data reduction program 'CxX4431'‘. 
It is listed as outlet air 
angle. Because of a difference 
in definition, it is necessary 
to use the outlet air angle 
Minus the stagger angle for 
BETA2. 


This is the ratio of specific 
heats. For the Subsonic Cas- 
cade Wind Tunnel the default 
value works well. 


This 1S an array of dimension 
10 of tolerances corresponding 
to MS and NOZES. Each grid 
solution will proceed until its 
TOLS value is satisfied. 


Inet 





MESH1 


MESHN 


LAMDAO 


CHORD 


Integer 


Integer 


Real 


Real 


Real 


None 


None 


None 


None 


This is an index in the 
arrays MS and NOZES of the 
first mesh to be generated 
and/or used for the flow so- 
lution. For geometry genera- 
tion runs, MESH1 selects one 
of the grids to be stored, 
provided NOZES(MESH1) > 0. 


This is an index in the 
arrays MS and NOZES of last 
grid to be calculated for 
geometry generation runs or 
used for the flow solution. 
For geometry generation runs, 
MESHN = MESH1L. Subsequent 
flow solution runs then solve 
the case for all grids listed 
in MS and NOZES between index 
MESHL and MESHN. For flow 
solution MESH] normally 
points to the coarsest mesh. 


To insure equal spacing of 
grid lines over the entire 
mesh, 


NOZES(I) - 1 
RIGS wise) © whieh o 


where MESH] S I S MESHN and 
Mr 2,035 et pas (RK = 2 2 t 
grid lines are doubled be- 
tween successive grids). 


Stagger angle of the blade 

row in degrees measured from 
the throughflow direction to 
the blade chord line. (Clock- 
wise is negative.) 


True chord of the blade, in 
the same units as the blade 
coordinates. 


Blade spacing in the same 
units as the blade 
coordinates. 


Lg 





The following parameters are required only for the 
geometry generation run. QSONIC gives the user a choice of 
two grid (mesh) generators. The Electrostatic Analog grid 
generator is applicable to any blade shape and any value of 
Stagger or turning. The interpolation scheme grid genera- 
tor works for most blades except those with sharp leading 
edges where the edge radius is less than .5% of chord. The 
interpolation scheme allows the user to concentrate grid 
lines in areas of high interest around the blade. With the 
Electrostatic Analog grid generator no concentration of grid 
lines is available. For the cascade configuration used in 
this study, the interpolation scheme provided gross errors 
in the flow solutions, so the Electrostatic Analog grid 


generator was used with good results. 


Name Type Default Description 
NED Integer None Total number of body definition 


coordinates that are input. These 
are the X,Y points that describe 
the profile of the blade, with the 
first point repeated as the last 
PONE. 


KN Integer None KN is used to indicate which grid 
generator is to be used. KN = 0 
will call the electrostatic analog 
grid generator. For interpolation 
scheme, KN = number of body points 
on upper surface from the minimum 
to maximum X points, inclusive. 


NED and KN are required for either grid generator. If 


the electrostatic analog grid generator is used, no other 
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parameters are used. 


If this grid generator is used the 


value from NOZES that 1S used must be odd. 


The following parameters are used only if the Interpo- 


lation Scheme is used to generate the mesh. 


Name 


RLE 


RIE 


THETL 


THe TT 


CAMPER 


STABAC 


Type 
Real 


Real 


Real 


Real 


Integer 


Real 


Default 


None 


None 


UToo9 


Description 


RLE is the leading edge radius 
of the blade, with units the 
Same as the coordinates de- 
fining the blade profile. 


RTE is the trailing edge radius 
of blade. 


THETL is the camber angle of 
the leading edge in degrees. 

It is measured from the aero- 
dynamic chord to the line tan- 
gent to the mean camber line at 
the leading edge; clockwise is 
positive. 


THETT is the camber angle at 
the trailing edge in degrees. 
It 1S measured from the chord 
line to the mean camber line at 
the trailing edge. Clockwise 
1s positive. 


For blades whose chordline lies 
outside the blade profile, such 
as the DCA blading discussed in 
this report, extra grid lines 
Surrounding the blade are 
needed to interpolate the blade 
position. The truncated value 
of MS( )/CAMPER is added to 

MS( ). The maximum allowed 
value of MS( ) + MS( )/CAMPER 
1s 30 grid lines. These are 
the grid lines that enclose the 
blade profile. 


STABAC is used only for test 
cases where no blade shape is 
to be input. Default value is 
usually adequate. 


ie 








—— an | ee 


CHOP 


SMOOTH 


LEON LY 


SLP1 
SaeP 2 
SLP3 
SLP 4 


The 


required 


Name 


MINF 


Real 


Logical 


Logical 


Real 
Real 
Real 
Real 


False 


False 


He be 


If leading edge or trailing 
edge radius is less than 2% of 
chord consult Ref. 9 for 
Claritecatiton. “Nexmally, 

0.9 < CHOP < 1.0. 


For automatic addition of more 
blade definition points in the 
region of the leading and trail- 
ing edges set SMOOTE = .TRUE. 
This should be done for all 
blades except cusps and wedges. 


LEONLY = .TRUE. for smoothing 
about the leading edge only; 
this is used if the trailing 
edge 1S a cusp or a wedge. 
SMOOTH must also be .TRUE. 


These parameters control the 
concenteation of grid lines, if 
desired. The default values 
worked well for the DCA blades 
reported herein. For con- 
trolling the amount and loca- 
tion of the concentration 
consult Ref. 9. 


following logical and case control parameters are 


only for flow solution runs. 


Type 
Real 


Default 


None 


Description 


This is the Mach number at the 
upstream boundary. This can be 
determined from the non- 
dimensional velocity X output 
from the cascade tunnel data 
reduction program 'Redd 5' and 
the relationship, 


iat ee. = = sahaetn oy elle 
5 ff — (es Y 





ds 
x? 2 
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TOLS (I) is the tolerance for 
MS (I) and NOZES(I). There are 
three forms of input permitted. 


A) -1.0 < TOLS(I) < 0.0: Cal- 
culations of the flow solu- 
tion will proceed until the 
relative circulation error 


N 
e a 
“CALC” TEXACT | ¢ |wons (1) 
EXACT 
3 


TOLS values between -10 


and lon are typical 
values for grids. This 
method of input is avpro- 
priate only for lifting 
(non-symmetric blades) 
cases. 


B) 0.0) < @Me@BSH I) < 1.0: Cal- 
culations of the flow solu- 
tion will proceed until the 
average relative chance in 
potential is less than the 
absolute value of TOLS(I). 


(=) < |TOLS (I) | 
> /AVE 


Typical values should be 
between LO= and il 


Gy 1.054226 Gr) +: Calecula- 
tions proceed until the 
number of iterations equals 
TOms(r). 


Regardless of TOLS input, the 
SOlution fer each grid will 
stop after 300 iterations 1f 
the TOLS criteria has not yet 
been made. All three forms 
were used for the cascade con- 
figuration reported herein 
with no discernible differences 
in results. 
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OVEREL 
UNDERL 
SUPREL 
NOWREL 
NOTYET 
TEGARD 
DAIMP 
oii 


Tt 


ALLOUT 


QUASTI3 


NSTRM 


Real 
Real 
Real 
integer 
Integer 
Real 
Real 
Real 


Integer 


Logical 


Logical 


Integer 


e 
oo uv 


OrFMNNNRFR EF 
° Oe 
NO © 


_~ 
Le) 


- FALSE. 


.PALom. 


The default values of these 
parameters are adequate for 
flow conditions in the subsonic 
cascade wind tunnel. 


Number of iterations between 
intermediate printouts of re- 
Siduals and Mach number. The 
information controlled by this 
parameter is of limited value 
in comparing with measured 
data, so a value greater than 
10 reduces the amount of com- 
Puser Prelntout. For the study 
reported herein 40 was used. 


To list the flow quantities at 
all grid points in the last 
mesh set ALLOUT = .TRUE.. Un- 
less a very coarse grid is 
used, the output resulting from 
ALLOUT = .TRUE. would be ex- 
tremely voluminous and of 
limited value. Until the cas- 
cade is configured so it is 
possible to take data from 
between the blades, ALLOUT 
should be .FALSE.. 


QUASI3 = .TRUE. to activate 
streamchannel thickness and/or 
radius variations. The cascade 
wind tunnel has no radius 
variations, but to simulate 3-D 
effects the streamchannel 
thickness is reduced at the 
exit boundary by a factor of 
1/AVDR. This data is placed in 
a file used by QSONIC if a 
quasi 3-D solution is desired. 


This is the position of desired 
Streamsurface data on the 
streamchannel file used if 
QUASI3 = .TRUE.. Currently the 
default value of 1 identifies 
the proper streamsurface data 
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RINF Real 
WAKE Real 
MINF2 Real 
OMEGA Real 
VAXIAL Real 
FLOCO Real 


die: 0 


Oi 
Oa 
@. 
Se 
oo. 


OOO CO © 


in the streamchannel file. If 
the output from the NASA code 
"Meridl' is used for the 
streamchannel data, then by 
using different values of 
NSTRM, different streamsurface 
data may be used. 


This is the spanwise radius 
at the upstream boundary di- 
vided by aerodynamic chord. 
Radius effects are activated 
if RINF # 1.0. The current 
version of QSONIC allows the 
following cases. 


QUASI RINF Results 
1) .FALSE. 1.0 Planar 2D flow 


2) .TRUr. #1.0 Thickness on 
file: radius 
on file. 


3) «TRUE. 1.0 Thickness on 
file; con- 
stant radius. 


Only 1) and 3) apply to the 
cascade wind tunnel. 


These parameters apply only if 
the test cascade is rotating 
and/or the downstream Mach is 
near 1.0. 


At this point all of the logical and case control para- 


meters necessary to use QSONIC for the flow conditions pos- 


Sible in the subsonic cascade wind tunnel have been discussed. 


The following is a description of the Bulk Data input for 


the geometry generation run and the flow solution run. 
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the format for all variables in the bulk data for mesh 
generation (geometry) is namelist form. The namelist is 


BNSTUPS . 


Name Type Default Description 





H2 Complex None This is a table of points de- 
fining the blade profile. The 
real vart = X, and imacinary 
part = Y coordinate. The table 
begins at the point of maximum X 
value at the trailing edge and 
proceeds clockwise back around 
to the first point, which is re- 
peatec. The blade must be at the 
stagger angle and the origin at 
the pGuint of Minimap XK. For 
Electrostatic Analog grids, the 
stagger angle must be positive 
(leading edge low, trailine edce 
high). The maximum number of 
points in H2 is 99 for the inter- 
nolation scheme or 63 for the 
Electrostatic Analog. 


BUG2 Logical .PRLoE. BUG2 = .TRUE. for a more detailed 
output of geometry generation. 
This will include the X,Y coordi- 
nates that define the mesh as 
well as second derivatives at 
grid points on the body. Except 
for trouble shooting this data is 
of limited value at the vresent 
time since there is no graphic 
Ole Pale. 


The bulk data for flow solutions consists of a mesh 
file and the streamchannel data file if quasi-3D effects are 
to be calculated. The mesh file is created by QSOMNIC during 


the mesh (grid) generation run. No further inputs are re- 


quired from the user for the mesh file. 
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The streamchannel data file must contain a table of 
streamtube thicknesses, radial positions and corresponding 


X values along the streamsurface. 


Name Type Default Description 





CHO Real None CHO is the aerodynamic chord 
multiplied by the cosine of 
LAMDAO. (LAMDAO = stagger angle) 


NRSP Integer None NRSP is the total number of data 
points in each of the tables of 
thickness, radial position, and X 
location. If MRSP is 2, a linear 
distribution is obtained between 
the endpoints given. NRSP = 2 
was used for the study reported 
herein with good results. 


RM Real Mone Array of corresponding X loca- 
tions for thickness and radius 
data values. xX = 0 represents 
the leading edge of blade, with 
the blade at stagger angle. The 
units can be any consistent 
length scale common to CHO, RM, 
RMSP and BESP. Inches were used 
in this study. 


RMSP Real None Spanwise radial positions of 
streamsurface at the X locations 
given in RM. RMSP was not used 
in the current study. 


BESP Real None Array of streamtube thickness 
values at the X locations speci- 
fied in RM. For the stucy re- 
ported herein, at X=0O0a 
streamtube thickness of 1.0 was 
arbitrarily selected. The 
streamtube contraction through 
the test section was similated 
by reducing the thickness at 
X = 0, by a factor of 1/AVDR at 
the trailing edge. Duval [Ref. 
3] explains AVDR. 





Dy. 3 


PrReeaRlNG INPUT FILES 


QSONIC was originally configured to use several input/ 


output devices while reading data and generating output. 


The input/output devices are listed below as used by OSONIC. 


P70 Unit 


2 


x3 


18 


23 


Usace 


File containing streamchannel data. This is used 
eniy WE QUASI = .TRUE.. 


Standard card input; NAMELISTS PARAMS, INSTUF. 
Standard printed output. 


For mesh generation runs, coordinates of all mesh 
points are written here. For flow solutions, X, 
Y, velocities, pressures and MACE are recorded 

for graphic display. The program currently has no 
graphic output capability. No user action is 
necessary to create this file. 


Used as temporary storage. No output is stored 
here. No user action iS necessary in conjunction 
“ath thys file. 


Previously developed mesh coordinates are read in 
from this file during the flow solution. After a 
mesh generation run, the user must create this 
Pee ama put in it the data from I/O unit 13, so 
that during the flow solution QSONIC can read in 
the mesh points. 


QOSONIC is presently configured to operate with the CMS 


system of the IBM 370 computer. This system provides a high 


degree of flexibility in parameter selection. With this 


system, all the input/output units previously mentioned are 


on the disk space assigned to the Turbopropulsion Laboratory 


eet) . 


Access to QSONIC and the TPL disk space can be ob- 


tained through the Director of TPL. 
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“hewtiust step in using QS@MIC is the creation of the 
data file necessary for the mesh generation run. This is 
done using the standard procedures of the XEDIT function of 
the CMS operating system. Reference 14 has specific in- 
structions for creating new files. The filename and file- 
type for the data used in this study was NAMELIST GFOMD. 
Once the data file is opened, the necessary data is input 
beginning at column 2 of the virtual card. Since the varia- 
bles are in namelist form, they can be input in any order. 

Table D.1 is an example of the data file necessary for 
the mesh generation run. The TITLE must appear on the first 
line (FORMAT = 20A4). After the TITLE, the logical and case 
control parameters are input after the namelist GPARAMS, 
When all the case control parameters required are input the 
PARAMS namelist is closed with &END. On the next line of 
the data file the bulk data for the mesh generation run is 
input in namelist form, with the namelist &INSTUF. The 
Hl = 100*(0.0,0.0) that appears after &INSTUF on Table D.l 
was used on earlier versions of QSONIC, but is not used in 
the present version. It should, however, appear in the data 
file before the H2 variables (X,Y coordinates defining the 
blade profile). 

At this point, some ciscussion of the coordinates de- 
fining the blade profile is warranted. Table D.2 is a 
listing of the X and Y coordinates of the DCA blading used 


im Gnas Stay. Figure D.1 is a plot of these coordinates. 
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Recall that the coordinates defining the blade profile for 
QSONIC must be for the blade at the stagger angle and mini- 
mam X at the origin. The coordinates of Table D.2 were 
translated and rotated using a coordinate transformation 
routine for the HP-67 programmable calculator. These new 
coordinates appear in the namelist INSTUFF on Table D.1l. 
Figure D.2 is a plot of the translated and rotated coordi- 
nates. It is highly recommended that such a plot be made 
for any new blade profiles, to ensure that the original 
coordinates are translated and rotated properly. 

The second step in using QSONIC is the creation of the 
Gata files necessary for the flow solution run. The file 
used for the flow solution in this report is on the "PL disk 
Space with a filename/filetype of NAMELIST FLOWD. The sim- 
plest way to open this data file is to use the YEDI™T func-~ 
tion, as discussed is Ref. 14, to start a new file. Then 
input the same data as is in the cata file for the mesh 
generation run usince the X™DIT subcommand GFT (filename) 
(filetype). The aporopriate changes anc additions can then 
be made to this file. Table D.3 is an example of the data 
file just discussed. 

Two more data files are recuired for the flow solution. 
The data for one of these is created by the mesh generation 
run. The other file contains the streamchannel thickness 


data for implementing quasi-3D effects. 





After the mesh generation run, a file with the file- 
name/filetype MESHGEN DATA will appear on the disk. Create 
a new file with the filename/filetyoe MESHIN DATA. This is 
most easily done by issuing the command 'XEDIT MESHIN DATA'; 
then use 'GET MESHGEN DATA'. This file contains the pre- 
viously developed mesh coordinates. 

The streamchannel data file should have the filename/ 
filetype of DATA5D DATA. The format for the data file is 


shown below. 


Virtual Card Column No. Variable Name 
1 BLANK 
2 BLANK 
3 21-30 GEO 
4 BLANK 
> 36-40 NRSP 
6 BLANK 
7 BLANK 
8 
9 
10 
Je | 
12 1=S0 (Sik. 5) RM 
¥ 


As needed 


1-80 (8F10.5) RMSP (not used in 


se this study) 


As needed 
1-80 (8F10.5) BESP 


Table D.4 is an example of the data file used in the 


present study. Since NRSP = 2 was used for this study a 
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lamear distribution is assumed for the streamtube thickness 
values and only 2 values of RM and BESP are required; there- 


fore, only 1 virtual card was required for each array. 


D.4 PROGRAM OUTPUT 

The output generated by QSONIC for the geometry genera- 
tion run includes a printed listing (I/O unit 6) and a mesh 
point file (I/O unit 13). The printed listing under the CMS 
system I/O unit 6 is normally the computer terminal unless 
the command ‘FILEDEF 06 PRINTER' has been invoked. It is 
unusual for the program to run properly the first time, so 
initially it is helpful to have the printed listing appear 
at the terminal. Once the program is running properly the 
output should be sent to the line printer. 

The flow solution run output consists of a printed 
listing (I/O unit 6) and a plot data save file (I/O unit 
13). Once the flow solution is running properly the printed 
listing should be sent to the line printer. 

Table D.5 is an example of the output generated by the 
geometry generation run. Figure D.3 is a plot of the grid 
output points on the blade surface, horizontal chord, pro- 
duced by the mesh generation run. Figure D.4 is a plot of 
the grid output points with the blade at the stagger angle. 
Table D.6 is an example of the output generated by the flow 


solution. 


PaO 





A detailed explanation of the printed output for the 


program QSONIC may be found in Ref. 9. 


D.5 RUNNING THE PROGRAM 
Mie fives on the TPL disk space that apply to QSONIC 

are listed below: 

OSONTC EXEC yal 

QSONIC FORTRAN Al 

Q@SOnNIC TEXT Al 

NAMELIST GEOM Pek 

NAMELIST FLOW Al 

NAMELIST GEOMD Al 

NAMELIST FLOWD Al 

DATAS DATA Al 

DATASD DATA Al 

QSONIC EXEC sets the input/output devices required to 

read and store data. QSONIC FORTRAN is the source program. 
To document the changes to QSONIC necessary to use the code 
with the IBM 370 operating system and serve as a reference 
for future users, a program listing is included at the end 
of this appendix. QSONIC TEXT is the computer executable 
code created when QSONIC FORTRAN is compiled. NAMELIST GEOM 
and NAMELIST FLOW are the data files for the geometry gen- 
eration and flow generation respectively for the example in 
Ref. 9. DATAS DATA is the streamchannel data required for 


the quasi-3D solution for the example in Ref. 9. 
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NAMELIST GEOMD is the data file for the geometry aqen- 
eration for the DCA blading used in the studv reported 
herein. NAMELIST FLOWD is the file for the flow solution 
for the cascade configuration used in this study. 

QSONIC exvects the input data to be ina file on the 
TPL disk space named NAMELIST DATA. Since the first time 
QSONIC is run is to develop the body centered mesh, the file 
NAMELIST GEOMD must be renamed NAMELIST DATA, using proce- 
dures specified in Ref. 13. Because QSOMIC requires large 
amounts of virtual memory, extra storage must be defined for 
the code to operate. This is accomolished by issuing the 
command 'DEFINE STORAGE 1504K'. 

With the data file renamed and more storace defined, 
type 'QSONIC' to load the program. The output will appear 
on the terminal screen unless FILEDEF 96 PRINTER was invoked 
pemer to loading the program. 

After the mesh generation is complete, rename NAMELIST 
DATA to NAMELIST GEOMD and change NAMELIST FLOWD to NAMELIST 
DATA. Create a data file with the filename/filetype MESHIN 
DATA. The elements of this file are the same as the elements 
in the file MESHGEMN DATA that was created by the mesh gen- 
eration run. The necessary input/output files are now con- 
figured for a flow solution run. Issue the command 'QSONIC' 
to begin execution. 

If the program output appears at the terminal it is 


possible to have some I/O error messages appear with the 
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output. This is because the write statements in QSONIC are 
formatted for the 132 character long line of the printer. 
These errors do not affect the validity of the program 
eutput. 

The explanation for any error or condition message 


generated by QSONIC can be found in Ref. 9. 


B-6 QSONIC UPDATE 

Recently an improved version of QSONIC was reported by 
NASA Lewis Research Center [{Ref. 13]. The new version re- 
quires less virtual memory and executes approximately 30% 
faster than the version presently in use at NPS. Also, the 
output appears in a different format than is described in 
this appendix. Reference 18 describes the most recent ver- 


sion of QSONIC in detail. 
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INPUT DATA FOR MESH GENERATION RUN OF QSONIC 
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TaAeee D.2. TEST BLADE COORDINATES 
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i Sa 06271 0.584 
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2 AG22 02272 0.671 
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INPUT DATA FOR FLOW SOLUTION RUN OF QSONIC 
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TABLE D.4. DATA FILE FOR QUASI-3D SOLUTION 


STREAMLINE DATA FOR DCAI1 
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SAMPLE OUTPUT FROM MESE GENERATION RUN 


TABLE D.5. 


uo a @¢ © 
JIN _ r) 
WATIPOo 
eNO ext 
eo z= G " i 
(v) —~— << © a OW 
® ad if © Wow 
uN HL © exo 
o PHS WO st 
gf LISO 2s2eo ew 
5 ef «© OO 
N ® «OS 
- © 
y © e 
aR) Hu. es NS 
za aH DAZ © 
Uy Va) Os © ome © 
© Pree HO 
OW) @ AS 
>» KO fe Ae) 
ew ao a ae) 
= ro OOO 
Uy Jr ly «© 
= wy © -O 
>) =Zemo Or O 
ud ee HO NN © 
— Ar Ooo Oo 
= HOO iL & 
= nw“ WoO = ZF e 
NY wHOoO mm © 
UJ Lcroe Co 
an) oO 6 6 » 
< Ox a 
ait ® iT) 
fe 8} ~ N 
" Uu. 
= ¢ us "4 ZZ 
b2 uoW Qa —_— 
Yt <i Ww = 
vs OM > aad - 
<< VINIW)D <I @ = 
SY) =Ouw - oOo 
<f{ -g=— 7a) a a 
= fier © ~ i 
ooafy OzO 
~ © e fo 
e © = & 
— OO MOS 
“uw& e zO 
© eS 
‘) iL 6 
2) Cl " O 
= a Own mo 
7S) oud rane OY 
UJ —«O NY 
Gee WwW fF We 
= NO ew = A - 
oD LJ tL i wcO 
Oa. = WAM Od eO 
=< UJ ‘an > ao 5 eFSiLO 
Tm fe LL <{ J CW aa 
° - ZOMN - 
ot J ui ee @&OH = ee 
G " © 9) © 
aw (od Oo uid 
Ore <— Oo as 
—— -— Oo = § 
MC Ye eONOM O = 
mW WOooCoeom dO 
uu Cc O00 e.« Oe 
> Ul ~ oOoCOax Cc 
VW -— OOo © 
Oat wil Oai©@ © 
= fs Oro il © 
za <I O Oo Oa © 
© x ad on Oe CO 
N «tL MPOUS ea e 
or | Qa©® <{ 
Ia (a ) 
=) 


GEND 


ELECTROSTATIC ANALOG GRID GENERATED 


20 


it 


49 IMAXY 


MOS ONS BOOP MOD OOP MO eudnor Qe wmoQ 
RAIN SAN FOO SHO OrHAIWDO OD OO MMM ODO 
ND et BOM Mh AMON OOD HO OH OOH Mea N a MOU 
TOR P= COPD ODP GND MNS tO DR MmNOUIST TP DOOU 
QM WOOP NODANSIKNAMNKHOMN NODOD Qed HO 
Og Or PSRO NDE D SOF NOP DD OOD OD MOUS 
ee a es Vg SOI NNN NAN INN 
are 0@ @ 6 ee 8etOs8ClUhthlU em CUD @e6e¢e8?@ e®*ees%eeee es @ 


ey lololololeleloiatelocoiererelor— OOOOOOCoeoe 


VS) 
WU 


OS MAA aT HID Smt MLA OO FOMMO OPO MAD 
LN BO AULD Pe Phe PDA MOD ah AHOOD Pe fO fh Pr SAO 
mt PP OND MLN OLN OR ON DD PENA Ut oh CY 
A BOF PAO AND ot FP OM OOSMONDAWOP> HAO 
CDN CO ed tt OO SB DONO OP UD ee ONION AO WO 
OD OOP mt LV ON MAS MO AF SHLD O UNO oF OD PO OU ~ 
LOOO OS 4 tt NN OM A PP PUN LAO OP PED OAT HAO 


wDeeteertreeeeeseeeseereee eevee e eee 6 @ @ 


ANDOOT DOD CODD CCQ COC COOOOOC OOo ooo 


NEGRI 


Uae: 





(Continued) 


TABLE D.5 


6 CHO 0.969 
URFACE, CHORD AT ST 


a0 
S 


LU 


Det YO NDP HOO ONO SND $F OO 
FLAP D SAQA OM ett PPOMmOLS 
MAA YVOSMAPPNYODOPANMINSMD SHO WO 
SVR POV a DS at SI POOP CO — 
NN St DIOPRDOOVWVMATHOOaAS ff a 
FPL AMAA HOD DOMOLPMHODOSL = 5 
MICIQUIQUNIOQIOS OI OI ad ed tS St Ht HOOOOOW 
°ee everest? Gee ee @€ &© & © © 8 eB 8 Of 
ek oo es SSO OO OSM ISLS 
Pomme 
‘enTo 4 
a WwW 
ONMDND = MAM OWSFOOWMDADOD fF 
DO = HOD OM DOD HHO AN DONS PFO Zz 
CO) FD POPU oe mt CP fe OO OOO at OOO 
PFOA tt DU OH OF DOD MS OM QOAWWOOD 
ODFPHDPALPOOMD DO Pea Aroayryananmac 
OWS OM OD ITOMN HP CID FOUONDOSTH © 
PO DMM OOQUAST FEO MANNN HOO OW LL 


oe ee eveveeeeeveet oe © © © © © @ ef) 


DOVBOOCOCDDOOOCOOCCOCOCOOOO = 


< 


== 


w59 


eleleleleleleleleleleleleleleleloer ff. | 
elelelelelelelel@lelelel@lelelelelelele) 
HPP HE HH HHH HHH HHH +H I SI 
CRS ULPULP OLS LAU A A 
RWOPDAOAAPANINODD LOND ON ovr 
CIN OD SD MAUI  HORIMOOO DOAN 
AID FOOD Fe DD ADD Re OF ed fe ee CU 
NAINA ated et St ted Sd HI OOD oP OO 
ed deed el peed ed eed ved ped med pend aed peed ead ed ed ed CY GY OO 
eeeeee%ee#sr#eege+@eg#sr#eseee#e##?eoe @©8 8® @ 


eleleloleleleleleleleleleleleolelolejlelo! 


elolelolelelelele\(elelelelolololelor |.) 
lolelelelelelelelslelelelelololelelelo) 
HHEHRERHE HEHEHE HH HH HHH I | 
UU UU OO OO OO OO dd 
OR SP CD SP COLD OAS Pe tO OO 
POM DO DONO OU F ONO MOS 
PDR TO POM SR AOS EMMA HHO 
DODD SF NODNMOM -eONN NA 
PPP FPS SS MMM MAUIAUN Sete OO 
*oe¢ee<c#ecfe30eo0Umr!ftmUmmlUC OWUCUMOWUlUCOWlCeklUCOlUCUCOMmUCOlhUC DMC‘ UCD hl HhlUS 


QDOOCOSCSO GOOG OOO COSOOGO0O000 


OM NOAPOOM DHROANMPIIOM™ OO 
wo td de 





TABLE D.5 (Cont‘d) 
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